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>- Abstract: We consider charged black holes in Einstein-Gauss-Bonnet Gravity with 

^ Lifshitz boundary conditions. We find that this class of models can reproduce the 

CO anomalous specific heat of condensed matter systems exhibiting non-Fermi-liquid 

behaviour at low temperatures. We find that the temperature dependence of the 
Sommerfeld ratio is sensitive to the choice of Gauss-Bonnet coupling parameter for a 
^ given value of the Lifshitz scaling parameter. We propose that this class of models is 

dual to a class of models of non-Fermi-liquid systems proposed by Castro-Neto et.al. 
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1. Introduction 

The AdS/CFT correspondence conjecture [1] has motivated the development of dual 
gravity models to describe strongly correlated systems in condensed matter physics 
(CMT) [2, 3, 4, 5]. This "AdS/CMT" correspondence, while more speculative than 
its AdS/CFT predecessor, has yielded some tantalizing insights into the behaviour 
of these systems. Systems of heavy fermions [6] - in which electrons near the Fermi 
surface of certain materials have an effective mass much larger than the free electron 
mass - are of particular interest since their behaviour is not that of a conventional 
Fermi liquid. The ratio of the specific heat to temperature (the Sommerfeld ratio 
denoted by 70), for low temperature, does not become constant but instead rises 
as the temperature decreases [7]. Indeed, there are some observations confirming 
that this ratio has logarithmic temperature dependence [8]. On the other hand, 
Castro Neto et al. proposed a model which predicts that the physical properties of 
the f-electron compounds at low temperatures are characterized by weak power law 
behavior, 70 = C/T oc T~^^^ where A < 1 is some characteristic parameter [9, 10]. 
This has been observed experimentally in a number of non-Fermi liquid materials 
[11, 12]. Though a full understanding of quantum criticality remains an outstanding 
problem, it is believed that such non Fermi liquid behaviour is due to the existence 
of a quantum critical point [13]. 
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Near a critical point many condensed matter systems are described by field 
theories with anisotropic scahng symmetry of the form 



t XH, 



X — )■ Ax, 



(1.1) 



where z{> 1) is a dynamical critical exponent represents the degree of anisotropy 
between space and time; manifestly z = 1 pertains to relativistic systems. Applying 
the AdS/CMT conjecture entails consideration of spacetimes whose metrics have the 
asymptotic form 



where the coordinates {t, r, x*) are dimensionless and the only length scale in the 
geometry is i. Such spacetimes are conjectured to be the gravitational duals to field 
theories with the scaling properties given in (1.1) [14]. It is straightforward to show 
that metric (1.2) has the scaling properties in (1.1) provided r — A^^r as well. 

Recently it was shown that these gravitational duals yield non Fermi liquid 
behavior for the specific heat of quantum critical system that is qualitatively similar 
to what is seen in experiments on heavy fermion systems [15]. The effects of finite 
temperature are induced by the presence of a black hole with the asymptotic form 
(1.2), and the addition of electromagnetic charge introduces another energy scale 
that can lead to interesting dynamics. The specific heat was found to scale as T"^/^ 
at high temperatures T, whereas at low T the behaviour of the specific heat was 
sensitive to the scaling exponent z. For z = 1 conventional Fermi liquid behaviour 
was recovered (with the Sommerfeld ratio becoming constant), whereas for z > 1 the 
Sommerfeld ratio increased with T as is observed for several heavy fermion systems. 

We consider in this paper the effects of adding in higher-order curvature correc- 
tions in the form of a Gauss-Bonnet term on Fermi liquid behaviour. In the context 
of AdS/CFT correspondence, inclusion of higher powers of the curvature is necessary 
in order to understand CFTs with different values for their central charges. In the 
context of condensed matter holography, they modify a number of important quan- 
tities, including viscosity, DC conductivity, the superconducting phase transition, 
and more [16, 17, 18, 19]. The simplest and most natural such term to include in 
the action is the Gauss-Bonnet term. Gauss-Bonnet gravity is a particular case of 
Lovelock gravity, which in (n + 1) dimensions yields p-order curvature terms (with 
p < [n/2]) in field equations with not more than 2 derivatives of any metric function; 
terms for which p > [n/2] are total derivatives in the action and so do not contribute 
to the field equations. To obtain a (3-|-l)-dimensional GMT, for which a (4-1-1) di- 
mensional gravity theory is the conjectured dual, it is therefore natural to include a 
Gauss-Bonnet term on the gravity side since for n = 4 we can only have p < 2. 

While we will ultimately be concerned with (3-|-l)-dimensional GMTs, we shall 




(1.2) 
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formulate our considerations in [n + 1) dimensions, making use of the action 

1 = j r^'x^ (^R-2A + aCcB - -/.uF'^" - \h^uH^'' - ^B^B'^^ , (1.3) 

where A is the cosmological constant and £gb = R^u-ysR'^'^^^ — 4:R^^R'^'^ + R^ is 
the second order Lovelock (Gauss-Bonnet) Lagrangian. The coefficient a has the 
dimension of (length)^ and can be positive or negative, F^^, = d[f^A^] with the 
Maxwell gauge field, and if^^/ = d[^Bu\ is the field strength of the Proca field 
with mass = C. 

Neutral black hole solutions with a = have been obtained both analytically 
and numerically [20, 21, 22]. The a = charged case in (n + 1) dimensions has been 
more recently considered [15, 23, 24]. In addition to obtaining a holographic descrip- 
tion of a strongly coupled quantum critical point in 5 dimensions with asymmetric 
scaling that exhibits the aforementioned anomalous specific heat behaviour at low 
temperature for heavy fermion metals [23], a more general study [24] of the d = 
version of (1.3) found a broad range of charged black hole solutions and a general 
thermodynamic relationship between energy, entropy, and chemical potential. So- 
lutions to neutral Lovelock-Lifshitz black holes were also recently obtained [25] and 
their thermodynamic behaviour for spatially fiat cases was found to be the same [26] 
as the p = 1 Einsteinian case [22]. Perturbative corrections to the charged Lifshitz 
case due to a Gauss-Bonnet term (a ^ 0) have been computed [28], but full black 
hole solutions and their implications have not been studied. 

We begin by introducing (n -|- l)-dimensional equations of motion in sec. (2), 
where we find the constraints that must be imposed on the constants in (1.3) in 
order that the asymptotic form (1.2) hold. We then investigate solutions at large 
r in sec. (3.1), which provide us with useful information about possible solutions 
due to different choices of a. Near horizon expansions are computed in sec. (3.2), 
and these are used to assist the numerical methods we employ to find solutions. We 
analyze the thermodynamic behaviour of these charged Gauss-Bonnet black holes 
numerically in sec. (4) where we specifically focus on z = 2 charged black holes in 
(4 -|- 1) dimensions. We find that the Gauss-Bonnet coefficient a plays a role dual 
to that of the characteristic parameter A in the Castro Neto et al. model [9, 10]. 
In fact, if the coefficient a is smaller than a certain minimal value am-, non- Fermi 
liquid behaviour is enhanced, and the temperature dependence of the Sommerfeld 
ratio 7o at low temperature can be characterized as a weak power law. However 
as a increases above the power law dependence of this ratio on the temperature 
becomes vanishingly small, similar to the A — )■ 1 limit in the Castro Neto et al. model, 
which as they claim can be fit with the logarithmic behaviour observed in some non 
fermi liquid compounds. The Gauss-Bonnet term restores Fermi liquid behaviour for 
a » am- We propose that Lifshitz-GB charged black holes are gravitational duals 
to the Castro Neto et al. model. Further numerical calculation reveals that the same 
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qualitative behaviour happens for arbitrary z > 1. We finish our paper by giving 
some closing remarks pertinent to our results. 



2. Field Equations in (n + l)-dimensions 

Using the variational principal the field equations that follow from the action (1.3) 
are: 

G^, + aG^^ + Ag^, = T^,, (2.1) 

V^i/^. = GB^, (2.2) 

di^B,] = H^,, (2.3) 

V^F^, = 0, (2.4) 



where 

T^u = -- {^Fp^pp'^g^y - PP^Fpy + -Hp^H'^'gpy - Hf^^Hp^ + G -BpE^gp^^ - BpB^ 

(2.5)- 

is the energy-momentum tensor of gauge fields, is the Einstein tensor, and G'^^ 
is given as: 

G'^^ = 2( — -R/^o-KT-R'"™!/ — '^R^puaR^'^ — '^RfiaR^u + RR^w) — 2'^GBfi'/ii/, (2.6) 

The (n + l)-dimensional metric that preserves the basic symmetries under consider- 
ation can be written as 

ds^ = t (^-r'^f\r)dt' + ^^^^ + (2.7) 

which is that of a black brane or (if appropriate identifications are carried out) a 
toroidal black hole. 

The gauge fields are assumed to be 

At = ir'K{r), Bt = qer'f{r)j{r), Htr = qizr'-^g{r)h{r)f{r), (2.8) 

with all other components either vanishing or being given by antisymmetrization. In 
order to get Lifshitz geometry (1.2) as the asymptotic form of the metric (2.7), one 
should demand /(r) = g{r) = h{r) = j{r) = 1 and k{r) = as r goes to infinity 
which impose the following constraints over the constants: 



_ in-l)z 2 _ 2(z - 1)L' 



& ' ^ z& ' 

[{z - 1)^ + n{z - 2) + n^] + n{n - l)a 

^cff = , (2.yj 
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with the following redefinition 
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-2a, a = (n-2)(n-3)d (2.10) 



where according to eq. (2.9) should be positive which implies a < £^/2 and 
consequently the cosmological constant Aog is always negative. Note that upon 
setting z = I'm (2.9) we recover an effective cosmological constant dependent on the 
coupling a rather than the usual AdS cosmological constant, i.e. A = —n{n — l)/2(?. 
Applying the ansatz (2.7) to the equation (2.4) yields the solution 

{r^^^y = ^Jg (2.11) 

where Q is an integration constant related to the Maxwell charge and we have chosen 
boundary conditions such that the Maxwell vector potential vanishes at the horizon. 

Substituting (2.9) and (2.11) into eqs. (2.1-2.3), the field equations reduce to 
the following system of first order differential equations: 

'■fr = 4(„-l)(fV-2.) '^ [(" - ^)(^ - ^)^-^ - - ^'"^ + - + "(^ - ^) + 

+ 2a{n - l)(n/ + n + Az - A) - 2{n - l){n + 2z - 2)e^g^ - QV^-^V}, (2.12) 

4r = 4(n-l)('v-2a) '^ K" ' ' + ' ^'"^ ' ' ' ' ' '''' 
- 2an{n - 1)(1 + g^) + 2n{n - l)fg^ + Q^^-^'^g^}, (2.13) 

4r = - 4(n-l)(k--2a) '^ K" " ^'<^ " ^'^^ " " ^'"^ ^ " + " ^' + 

+ 2a(n - 1)(V + n-A) -2{n-l){n- 2)fg^ - Q^r'^^^^'g^} + zgh, (2.14) 

= (n-l)(j(7-/i). (2.15) 

3. Series Solutions 
3.1 Solutions at large r 

Since we wish to understand the effects of the Gauss-Bonnet term on the previously 
considered [24] (n -|- l)-dimensional charged Einstein solutions with arbitrary z, we 
begin by linearizing the system in (n + 1) -dimensions. Requiring the general metric 
(2.7) to asymptotically approach the Lifshitz one, we write 

/(r) = l + wfi{r), 
g{r) = 1 +wgi{r), 
j(r) = 1 + wji{r), 

h{r) = 1 + whi{r) (3.1) 
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where we note that it is also necessary [20, 21] for j(r) and h{r) to approach unity 
in order to obtain (1.2). 

In previous section, we aheady introduced equations of motion which are shown 
in (2.12-2.15). Since /(r) does not contribute to the equations for g{r), h{r) and j(r) 
let's first study a set involving {(?, If one inserts perturbative expansion (3.1) 
into equations (2.13-2.15) one obtains the equations for small perturbations as 

n —z{z — l)/(n — 1) 

1 — n n — 1 

M3i{n,z,a) -z{n + z -2)/{n - 1) 

+ 4(n -l)r-'-^ (") 
which is identical to the Einsteinian case except for the term 

where we have rescaled r r/i. We have also included the Maxwell gauge field as a 
first order perturbation, which means that we substitute Q^/r^*^""^) = wQ'^ / r'^^"'~^\ 
since its falloff may be slower than other terms in the metric functions. 
Solving the set (3.2) one obtains: 

( ) = ^1^1 _ CsCh {n - 2z)eQ' 

j.z+n-1 ^(2+n-l+V7)/2 ^(z+n-l-^)/2 AA{n — Z — l)r'^"'~'^ 



hi{r) 
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^(2+n-l+^)/2 ^{z+n-l-^)/2 

[(2n-z-2)e-2{2n-3)a]eQ' 



■ / \ _ CiJi C2J2 CsJs 

Jlv) — ~ ^z+n-l ~ ^(z+n-l+^)/2 ~ ^(z+n-l-^)/2 

[{n + z- 2)i^ - 2(n + 2z- 3)a]i^Q^ 



(3.6) 



where 



((n-3z)'^ + 6n-2z-7)i^-2(n^-6nz + 6n + Uz + z^-15)a , , 
7 = -j2 ' (3-7) 

A = (n + ^ - 2)(n - ^ - \)e - 2{n^ - 3n - ^ + 3)a (3.8) 

C = z{z-l)L^ =^ r =— 

' {n - - 2{n - l){n + z - 2)a' ' n-V ' n-V ^ ^ 
^ z{n + z- 2)e - 2{n + 2z~ 3)a 
' (n-l)2£2_2(n-l)(n-z + 2)a' ^ ' 

J n-3z + l-^ n-3z + l + ^ 
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Substituting above expressions into the equation for small perturbations of /(r) 

rf^,, {n-2 + 2z)e^ -2{n + Az-A)a z{z - 1) ^ . . 

r^^hir) = g,{r) - -^^h,{r) 

+ (^-1)^1(0- ,^. (3.12) 



we obtain 
with 



CiFi C2F2 



[^n + z-2)e-2{n + 2z-2>)a]eQ^ 

4L2A(n - l)r2"-2 ' ^ ■ J 



_{n + 2z- 2)f - 2{n + 4z - A)a z-l 

- {n + z-l)L^ + n + z-l' 

z(z-l) 



{n- l){n + z- 1)' 



(3.14) 



{n + 2z - 2)e^ - 2{n + Az - A)a 2(^-1) 
'~ in + z-l + ^)L2 + + ^ ' 

2z{z - 1) 



(n - l)(n + z- 1 + ^)' 



(3.15) 



^ _ (n + 2z - 2)£2 - 2{n + Az - 4)a ^ ^ 2(z - 1) ^ 
^ (n + 2; - 1 - -^1^)1^ ^ n + z — 1 — ^ ^ 

2z{z - 1) 



{n-l){n + z-l-^y 



(3.16) 



There are three eigenmodes, of which two are decaying. The third one might be 
decaying, growing or independent of r depending on the value of a. There is a zero 
mode if 

«(o) ^ - ^ - 1)^' (3 17) 

2(n-2) ^ ^ 

Accordingly, for a. > a.'^^^ in order to approach the Lifshitz fixed point asymptotically, 
fine-tuning of initial values is required when numerically solving the non-linear field 
equations. However if a < q;^°\ all three eigenmodes are decaying, yielding a family 
of solutions with the same event horizon that asymptotically approach the Lifshitz 
background with different fall-off rates. 

3.2 Near horizon expansion 

To investigate the near horizon behavior of the solutions we consider the following 
expansions 

f{r) = /oV^^(l + Mr - ro) + /2(r - ro)' + •••), 
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9{r) = I ^° (1 ^g\{r - To) + g2{r - Tq)^ H ), 

i{r) = joVr - ro(l + ji(r - tq) + j2(r - tq)^ H ), 

/i(r) = /io(l + /ii(r - ro) + /i2(r - ro)^ + ■■■), (3.18) 

where tq is assumed to be the horizon radius. Inserting this into Eqs. (2.12-2.15) and 
demanding the coefficients for each power of (r — ro) vanish yields relations between 
the various constant coefficients in terms of ro and ho. For example go is given by: 



_ . 1/2 / 2(72 - 1) 



where 



gL = + - 2)n +{z- If) e - [n' + (2. - 3)n + 2(. - 1)^) i^'m) 



is the upper bound for the Maxwell charge for a given a. Since a < P jl the right 
hand side of eq. (3.20) is always positive. Clearly Qext decreases as a increases. 
Furthermore, from (3.19) we see that real series solutions exist provided 



with Q < Qcxt- AH other constants in the series solutions (3.18) can be obtained in 
terms of ro and ho as well but we don't write them here because they are too lengthy. 



4. Numerical analysis and Thermodynamic Behaviour 

To find numerical solutions for the system of ODE's (2.12-2.15) we apply the shooting 
method, adjusting initial values for the fields /(r), g{r), j{r) and h{r) and requiring 
each to approach their asymptotic values of unity. To find these initial values we 
use the near horizon expansions (3.18) by substituting ro -|- e for r where e and 
choose some values for ro and ho such that the fields approach unity at large r to 
within a certain tolerance of 10~®. While we can obtain explicit numerical solutions 
for the metric functions, we are interested here in their thermodynamic behaviour, 
particularly insofar as how the Gauss-Bonnet term modifies the Sommerfeld relation. 
Throughout we restrict our numerical analysis to 5 dimensions. 

4.1 Thermal behaviour 

The Hawking temperature of a general black brane (2.7) is 

T=^r,(ro,g,«), (4.1) 
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Figure 1 I T2 vs o. in 5-d.im6iisioiis with, z — 2. Loft; riGutral black holes with (yrn 

^ 0.22. 

Right: weakly charged black holes {Q 0.676 Qc) with am ~ 0.3. We have set £ = 1. 
with 

= ^, zy^l (4.2) 
9o 

where /o and Qq are (rg, Q, a)-dependant coefficients from the near horizon expan- 
sions (3.18). Noting from eq. (2.12) that there is rescahng symmetry for /(r), we 
write /o — t- /o/y^- Defining Q = Q/r^"^ we find numerically that is just a 
function of {Q, a) and does not depend on tq explicitly. 

To examine the thermal behaviour of these charged Gauss-Bonnet black branes 
we follow two different approaches. First we investigate how temperature depends 
on the Gauss-Bonnet coefficient a for fixed electric charge. Then we investigate the 
variation of temperature with electric charge for fixed Gauss-Bonnet coefficient. 

Consider ffist fixed electric charge Q. For z = 1 the temperature function ti 
depends linearly on a as 

while numerical analysis indicates that the situation is different for z > 1. The reality 
condition from eq. (3.19) imposes the constraint a < a^xt, where 

_ 2n' + 2{z-2)n + 2{z-iy-Q' , 

2n^ + 2{2z-3)n + 4{z-iy ' ^ ' 

RecaUing that a < i'^/2, we define Qc = Q(aext = ^^/2) from (4.4), yielding 

Ql = n{n-1). (4.5) 

For charged black holes with Q < we have a^xt > ^^/2 implying the bound 
a < f/2; otherwise, if Qc < Q < Qcxt(= <5cxt/ro"^) then a < Ocxt < ^^/2. 
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Figure 2: Functions describing 5-dimensional neutral Einstein (Q = 0, a = 0) in red 
lines and neutral Gauss-Bonnet (Q = 0, a = 0.36) in blue lines both with z = 2. Left: 
metric functions /(r)(solid) and <^(r) (dashed). Right: gauge functions j(?')(solid) and 
/i(r) (dashed). We have set £ = 1. 

Numerically solving the field equations, we find that weakly charged {Q < Qc) 
or neutral black holes with z > 1 get colder as a decreases from its maximal value 
of £^/2 down to a specific value, denoted am, which can be determined numerically. 
As a further decreases below the temperature starts increasing. In other words 
the temperature is minimized at a = < ^^/2. 

We depict this behaviour in fig. 1 for small and large black holes in 5 dimensions 
with z = 2. It is evident from fig. 1 that it is possible to find isothermal pairs of 
black holes. For example, a z = 2 neutral Einstein black hole {a = 0) with tq = 0.4 
is isothermal to a neutral GB black hole with a = 0.36 of the same z and size. We 
illustrate in fig. 2 that these are indeed two different solutions. 

For charged black holes with intermediate values of the electric charge - Q ^ Qc 
- we find that the temperature has a minimum nonzero value at a < £^/2 and then 
increases monotonically as a decreases from i'^/2. There are no isothermal pairs. We 
illustrate this behaviour for typical cases in the left-hand side of fig. 3. 

The temperature T of strongly charged GB black holes, where Q > Qc (and 
z > 1) exhibits markedly different behaviour. In this case a < acxt < ^^/2 and we 
find, as shown in the right-hand side of fig. 3, that T starts from zero at a = Oext and 
then monotonically increases with decreasing a. Again, no isothermal pairs occur. 

Turning next to the question of how temperature changes with charge for fixed 
a, we shall consider separately the high temperature and low temperature behaviour. 

For a given a, high temperature behaviour corresponds to Q — >■ 0. In this limit 
the temperature function Tz is just a function of a and independent of tq. Hence 
T oc Tq for charged GB black holes at high temperature with arbitrary z. 

We have already noted that in the high temperature limit there exist pairs of 
isothermal black holes with different GB coefficients, one with a > am and the other 
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Figure 3: T2 vs a in 5-dimensions with z = 2. Left: intermediate charged black holes 
(Q ^ 0.902 Qc) with Oext ~ 0.55. Right: strongly charged black holes {Q ^ 1.128 Qc) with 
Oext ~ 0.425. We have set ^ = 1. 




Figure 4: r2 vs Q for z = 2 black holes in 5 dimensions with Om ~ 0.22. In the right 
hand plot solids show a > Om where dashed show a < am- 

with a < am- (The minimum GB coefficient am ~ 0.22, 0.02 and —0.1 for z = 2, 3 
and 4 respectively). We consider two such isothermal cases in fig. 4. We see that for 
small Q (or high temperature) the plots for the isothermal pairs a = (0, 0.36) and 
a = (—0.2, 0.4) (where am — 0.22) are essentially the same. 

However the distinction between pairs grows with increasing Q. As the low tem- 
perature limit {Q — Qcxt) is approached, the behaviour of becomes significantly 
different between the elements of the isothermal pair. For a < am the low tem- 
perature limit (T — )■ 0) is non-linear as a function of increasing Q, whereas for its 
isothermal counterpart with a > am, it tends toward linearity. We will revisit this 
distinction in the next section where we will discuss the behaviour of the specific 
heat at low temperature. 

Fig. 4 illustrates this for 2; = 2. The same situation occurs for arbitrary values 
of 2 > 1, as shown in fig. 5 for z = 3. We conclude that this is a general feature of 
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Figure 5: Left: vs a for 5-dimensional z = 3 neutral black holes. Right: T3 vs Q for 
5-dimensional z = 3 GB black holes with a = 0.18 (blue) and a = —0.24 (red) [am ~ 0.02]. 

planar Gauss-Bonnet black holes of arbitrary size and z > 1. 

In general for a given z, black holes with a < am show non-linear behaviour at 
low temperature. As a decreases below am, the maximum temperature, T{Q = 0), 
increases as does Qext- As a result the non-linearity at low temperature becomes 
more significant as a jam becomes smaller. However as a increases above am, the 
maximum black hole temperature increases while Qext decreases. In other words the 
temperature falls faster, with the low temperature behaviour becoming increasingly 
linear as shown in fig. (4-right). 

We close this section by noting that for a = am no isothermal pairs exist. 

4.2 Specific heat and Fermi liquid behaviour 

The specific heat at fixed volume is given by: 



is the Bekenstein-Hawking entropy of the black brane. A characteristic behaviour 
of a Fermi liquid is that its specific heat depends linearly on temperature. A recent 
study of Fermi liquid holography compared asymptotically AdS charged Einstein 
black holes(z = 1) to asymptotically Lifshitz ones {z > 1) [23]. The former obeyed 
Fermi liquid behaviour at low temperature, whereas the latter did not. Instead 
asymptotically Lifshitz Einstein-Maxwell black holes (2; > 1) at low temperature 
deviate from Fermi liquid behaviour. Indeed, for this class of black holes C/T diverges 
in the T — t- limit. 




(4.6) 



where 



S 




(4.7) 
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Figure 6: logiQ{C/T) vs log^gT in 5 dimensions with fixed charge Q = 1, z = 2, and 
am ~ 0.22. 



To see how the GB correction modifies specific heat, we begin with the z = 1 
case. Here we can easily obtain an analytic expression for the specific heat 



C 
T 



A{n - I)^7r2+V| 



r(f) (2n-3)Q2 + 2n(r2-l)(l-f, 



(4i 



using eqs. (4.1), (4.3) and (4.7). From eq. (4.8) it is straightforward to show that 
at high temperature (i.e. Q — > 0) the Sommerfeld ratio C/T oc T"~^. However at 
low temperature, where the electric charge approaches the extremal value Qcxt(= 
Qcxt/rQ~^) introduced in (3.20), the Sommerfeld ratio becomes constant: 



C 



TT2 



' ext 



T nVim- 



as 



T -> 



with 



ext 



(4.9) 



(4.10) 



_2n(n-l)(l-f; 

The low-temperature constancy of the Sommerfeld ratio is a key characteristic of 
Fermi liquids. 

We have already observed that T oc for arbitrary values of z in the high 
temperature {Q — )■ 0) limit. Consequently, since specific heat is proportional to 
r^-\ we find that C oc T^^-^)/^ for all GB black holes. 

The low temperature limit C/T for 2; > 1 must be obtained numerically. We 
have numerically computed the specific heat (4.6) for a range of values of a. In fig. 
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6 we plot the Sommerfeld ratio C/T as a function of temperature on a log-log scale 
for 2 = 2 for the isothermal pairs in fig. 4. We see although both a < a^n and 
a > am are isotherm pairs at high temperature, dramatic growth appears in C/T as 
T — )■ for a < am, increasingly so as a becomes much smaller than am- Therefore 
the low temperature behaviour becomes increasingly non-linear for GB black holes 
with a < am, with the Sommerfeld ratio depending on temperature as a weak power 
law, consistent with the prediction of the Castro Neto et al. model. The slope of the 
line in the log-log plot of Cy/T at low temperature from the bulk depends on the 
value of a, whereas for the Castro Neto et al. model it is characterized by A — 1. For 
example, for a = —0.2 the slope is approximately —0.43 (corresponding to A = 0.57) 
while it is approximately —0.38 for a = (equivalent to A = 0.62). 

The a > am isotherm counterparts exhibit rather different behaviour at low 
temperature, with the T-dependence of the Sommerfeld ratio showing vanishingly 
small power law behaviour as a increases (corresponding to A — )■ 1). For a » am, 
this becomes almost constant, restoring Fermi liquid behaviour. 



5. Closing remarks 

Adding a Gauss-Bonnet term to the Einstein-Hilbert action results in interesting 
low temperature behaviour for asymptotically Lifshitz charged black holes. By in- 
vestigating the thermodynamic properties of such black holes at low temperature, 
we found that the Sommerfeld ratio 70 is temperature dependent, governed by a 
weak power law behaviour depending on the value of a. This kind of behaviour is a 
characteristic property of some heavy fermion metals such as f-electron compounds 
at low temperature. A CM model of these non- Fermi liquids [9, 10] predicts weak 
power law behaviour for other properties such as heat capacity and susceptibility, 
i.e. 7o(T) = C{T)/T oc x(T) oc T-i+^ with A < 1 . 

We propose that Lifshitz-GB charged black holes are gravitational duals to the 
Castro Neto et al. model [9, 10], with the Gauss-Bonnet coefficient a holographically 
dual to the characteristic parameter A. To more firmly establish this duality will 
involve a study of other thermodynamic quantities like resistivity or susceptibility 
using the gravity theory. We leave this topic for future work. 
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